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Abstract 

We study a racetrack model in the presence of the leading a'-correction in flux 
compactification in Type IIB string theory, for the purpose of getting conceivable 
de-Sitter vacua in the large compactified volume approximation. Unlike the Kahler 
Uplift model studied previously, the a '-correction is more controllable for the meta- 
stable de-Sitter vacua in the racetrack case since the constraint on the compactified 
volume size is very much relaxed. We find that the vacuum energy density A for 
de-Sitter vacua approaches zero exponentially as the volume grows. We also analyze 
properties of the probability distribution of A in this class of models. As in other 
cases studied earlier, the probability distribution again peaks sharply at A = 0. We 
also study the Racetrack Kahler Uplift model in the Swiss-Cheese type model. 
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1 Introduction 



Recent cosmological data strongly suggests that our universe has a vanishingly small positive 
cosmological constant as the dark energy, 

A ~ W~ 122 M P (1.1) 

where Mp is the Planck mass []}j4] (and references therein). Arguably, its smallness is one of 
the biggest puzzles in modern fundamental physics. By now, it looks hopeless to search for a 
natural reason for its smallness within four-dimensional quantum field theory. 

A typical flux compactification in string theory involves many moduli and three-form field 
strengths with quantized fluxes (see the review pi). Together with the quantized fluxes, the 
moduli and their dynamics describe the string theory landscape. Since string theory has many 
possible vacuum solutions (i.e., leading to the so called cosmic landscape), it is argued that the 
spacing 5A can be exponentially small and possible values of A can have large ranges. As a 



result, such a small A (1.1) can easily be that for one of the solutions [6]. However, this alone 
does not explain why nature picks such a very small A, instead of a value closer to the string or 
Planck scale. 

Recently we proposed that a combination of string theory dynamics together and some basic 
probability theory may provide a natural order-of-magnitude explanation to this naturalness 
puzzle [7[j9]. The basic idea is quite simple. Suppose we can determine the four-dimensional low 
energy supergravity effective potential V for the vacua coming from some flux compactification 
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in string theory. To be specific, let us consider only 3-form field strengths F$ wrapping the three- 
cycles in a Calabi-Yau like manifold. (Note that these are dual to the four-form field strengths 
in 4 dimensional space-time.) So we have V(F^, <pj) — > V(rii, <pj), {i = 1, 2, N, j = 1, 2, K) 
where the flux quantization property of the 3-form field strengths F%s allow us to rewrite V 
as a function of the quantized values n« of the fluxes present and <pj are the complex moduli 
describing the size and shape of the compactified manifold as well as the coupling. There 
are barriers between different sets of flux values. For example, there is a (finite height) barrier 
between ri\ and ri\ — 1, where tunneling between V(ni, n 2 , n^, <f)j) and V{n\ — 1, n 2 , n^, 4>j) 
may be achieved by brane-flux annihilation |6]. 

For a given set of rii, we can solve V(ni,<pj) for its meta-stable (classically stable) vacuum 
solutions via finding the values 0j,min(^i) at each solution and determine its vacuum energy 
A = A(rii). Collecting all such solutions, we can next find the probability distribution P(A) of 
A of these meta-stable solutions as we sweep through all the flux numbers n^. Since a typical 
rti can take a large range of integer values, we may simply treat each rii as a random variable hi 
with some uniform distribution Pi(bi) and find the properties of -P(A). 

Let us focus on Type IIB string theories with known de Sitter (dS) vacuum solutions. It turns 
out that the string theory dynamics (i.e., the resulting functional form of A(6j)) together with 
simple probability theory typically yields a -P(A) that peaks (i.e., diverges) at A = 0. In fact, this 
peaking at A = behavior is relatively insensitive to the details of the smooth distributions Pi(bi) 
and becomes more divergent as the number of moduli/fluxes increases. However, this divergence 
is always mild enough so -P(A) can be properly normalized, i.e., J P(A)dA = 1, henceforth 
implying that the probability at exactly A = will remain exactly zero. Since the number of 
moduli as well as cycles that fluxes can wrap over in a typical known flux compactification is of 
order (9(100), a vanishingly small but non-zero A appears to be statistically preferred. In fact, 
it is not hard to find the likely value of A at comparable magnitudes as the observed value (1.1 ). 

Although the overall emerging picture is encouraging, there are many open questions that 
need to be more fully addressed. In this paper, we like to study an important issue in (7l|8l. In 



the Kahler Uplift model with the leading order a '-correction 10 13 in models similar to the 
Large Volume Scenario 14 , the large compactification volume approximation assumed is only 
moderately satisfied a posteriori for the SUSY breaking meta-stable solutions around A = 0. 
Since the large volume approximation works well in the presence of this and other a'-corrections 
as well as the stringy loop corrections 15 -20 , the constraint on the volume size leads to concerns 
on the validity of the approximation. 

The Kahler Uplift model studied has a single non-perturbative term in the superpotential 
W. To relax the constraint on the volume size, we generalize the model to include two non- 
perturbative terms in W, i.e., the racetrack model. (We like to point out that this has been 



briefly studied in 13,211.) Owing to this racetrack property, we find that the model admits 
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solutions with a large adjustable volume. 

Interestingly, in this Racetrack Kahler Uplift model, the stability condition for both the real 
and imaginary sectors requires that the minima of the potential V always exist for A > at 
large volumes. Further, the cosmological constant A is naturally exponentially suppressed as a 
function of the volume size, and the resultant probability distribution -P(A) for A gets a sharply 
peaked behavior toward A — > 0, which can be highly diverging. This peaked behavior of -P(A) 
is much sharper than that of the previous Kahler Uplift model with a single non-perturbative 
term studied in (7j|8). 

The paper proceeds as follows : The racetrack model is introduced and reviewed in section 
[2] Among the possible solutions for dS vacua, we focus on the set that allows large volumes 
that are not bounded from above. We contrast this with the single non-perturbative term model 
studied earlier. Although the range of A may be unbounded from above, these dS vacua are 
forced to have an exponentially small A in the large volume limit. In section [3} we present the 
probability distribution -P(A) for these dS vacua, which peaks sharply at A = 0. In section 
[4j we extend the racetrack potential to the Swiss-Cheese type model. Section [5] contains the 
discussions and remarks. Some details are relegated to appendix |A} 

2 A racetrack model and property 
2.1 Background 

We have focused so far on Type IIB models. It is important to comment on the difference 
between type IIB and IIA models with respect to the moduli stabilization. The moduli are 
four-dimensional light scalar fields parametrizing the geometric size and shape (deformation) 
of the compact six-dimensional internal spaces (as well as the dilaton-axion mode) in string 
theory needed to describe our effectively four-dimensional universe. The moduli stabilization in 
type IIA is typically very difficult to achieve since we have to stabilize the entire set of moduli 
simultaneously due to the absence of hierarchical structures. If we have no specific structure 
in the potential, we may expect that the mass (squared) matrix is given rather randomly at 
dS extremal points. Then the probability that all eigenvalues of the random mass matrix are 
semi-positive (required for meta-stability) is described by a Gaussian suppressed function of the 
number of moduli [22j[23] (see also [24]). Since we may expect typically C(100) of moduli fields, 
it is clear why type IIA stabilization is so difficult to find. On the other hand, we have the no- 
scale structure in type IIB; so the Kahler sector can be considered separately from the complex 
structure and the dilaton sector, which are stabilized at higher scales. As a result, the number 
of moduli to be simultaneously stabilized is drastically reduced. This hierarchical structure 
holds even if we introduce non-perturbative terms and a'-corrections accordingly, which weakly 
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break the no-scale structure so that a non-trivial potential is generated. Some models have few 
Kahler moduli and large number of complex structure moduli. Recently it is suggested that the 
large number of complex structure moduli helps to enhance the hierarchical structure [9]. So 
the IIB models are well-motivated to achieve the moduli stabilization with positive cosmological 
constant. 

As we just pointed out, a corner of the string theory landscape where moduli stabilization 
can be addressed explicitly is type IIB compactified on orientifolded Calabi-Yau three-folds. 
The four-dimensional effective action of the geometric moduli is given by a M = 1 supergravity 
theory of a set of chiral multiplets consisting of the dilaton-axion S, h 1,1 number of Kahler 
moduli Ti, and h 2 ' 1 number of complex structure moduli U. The past decade saw some progress 
for the stabilization of S and U from the use of quantized fluxes of three form field strength of 
the Ramon d-Ramond (RR) type F 3 = dC% and Neveu-Schwarz (NS-NS) type if 3 = dB 2 , which 
form a SL(2,Z)s covariant three-form: G 3 = F 3 — iSH 3 . They wrap three-cycles inside the 
manifold. Their 10-dimensional duals are 7-form fields wrapping dual three-cycles, which result 
in effective four-form constant quantized field strengths in four-dimensional space-time. The 
flux stabilization procedure operates supersymmetrically at a high scale. The Kahler moduli 
are typically stabilized at a parametrically lower scale than the complex structure moduli via 
perturbative or non-perturbative interactions. We set Mp = 1 throughout this paper. 

To be specific, we consider a racetrack model defined by the Kahler potential K and the 
superpotential W, 

K = - 21n \V + - In (S + S) - In (-i J O A ttj , 

V =(T + T fl\ i = - - £ 3 o ) — X (M) (S + S) 3/2 ~ -8.57 x l(T 4 x (S + sf 2 , 



71 



W =W (U i , S) + W NP , W (U h S)=Y, Jg 3 AQ = C^Ui) + iSC 2 (Ui) 
V =e K {K rj DjWDjW - 3 \wf} 



cycles ' 



where £ is the term coming from 0(a ) corrections to SUGRA 25 , and shows up as an uplifting 



term in the potential 10 u 12 . The non-perturbative term Wnp will be specified below, it is 
expected to be small compared to the tree-level flux contribution. Here the dimensionless 
compactification volume V = vol/ a' 3 is measured in string units. We see that the holomorphic 
3- form Q depends on the complex structure moduli C/j. 

A simplified version of Wo(Ui, S) has been discussed in earlier works [8], where one finds 
the supersymmetric Wq and then the ratio of the median value and the average value of the 
magnitude of A, |A|/ (|A|), which tends to decrease exponentially as the number of Ui increases. 
Since the stabilization of Ui and S are assumed to take place at a higher scale than that of the 
Kahler moduli, and this part of the analysis is very similar to the earlier work, we shall simply 
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assume a value for W$ in this paper and focus on the dynamics of the Kahler moduli. 



2.2 The racetrack model 



The non-perturbative terms in the superpotential W as in [26] (see also 27 29]) are crucial in 
the Kahler moduli stabilization. Compared to our earlier work [7j[8], the main new feature here 
is the presence of the new term Be~ bT in W . Together with the other terms in W, this forms 
the so called "racetrack". To focus on this feature, let us assume that dilaton and complex 
structure moduli are stabilized supersymmetrically D$W = dsW + {dgK)W = 0, D^W = at 
some higher energy scale, so W is reduced to 

W = W + W NP = W + Ae~ aT + Be~ bT (2.2) 

where the coefficients a and b are taken to be positive real. The potential V for the Kahler 
modulus T becomes 



V=e K 



K TT (d T WdfW + d T W{Wd T K) + d f W{Wd T K)) + ^ £ + J^ + ^l 9 \W\ 2 



(2.3) 

In this paper, we are mainly interested in meta-stable dS vacua achieved with the uplifting 
a'-correction term £ in the Kahler potential. 

Note that a solution of this type was found and used for an inflationary universe analysis 
in 21 . Here we are interested in more detail properties and systematic understandings of this 
class of solutions, including allowed region of parameters for dS uplift, such that we can analyze 
also the property of distribution of vacua. We also like to point out that the model has been 
briefly studied in [13], but in a different parameter region. 

Before proceeding, it would be interesting to see what happens for supersymmetric vacua 
obtained before turning on a'-correction. If the supersymmetric condition DjW\^ = o = holds 
before turning on a'- correct ion, the potential becomes 

V = e K (\D T W\ 2 - 3\W\ 2 ) ~ 3e K ^°\W\ 2 (-1 + |j (2.4) 

where in the last equation, we kept up to the leading order of £/V. We see immediately that dS 
vacua V\ m i n > by the uplifting from SUSY AdS requires £/V > 0(1). This clearly violates 
our assumption that a'-correction is under control in type IIB SUGRA approximation. Thus 
we need SUSY breaking vacua before the uplift by the leading order a'-correction. Note that 
the leading term of IDtW^ 2 starts with 0(£ 2 ) due to the supersymmetric condition^] 

1 Although our interest throughout this paper is for the uplift by the leading order a'-correction, it is worth 
commenting that the SUSY vacua can be uplifted to dS by introducing an explicit SUSY breaking term like 



D3-D3 pairs contribution considered in 26 30 . 
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In the large volume situation where (/V < 1, which is of interest here, the potential may be 
approximated to 

x ( x ^) = -^ cos y--— cos (f 3 y) + ^9j2^ 



X* Z X" X 

x = at, y = ar, z = A/B, = b/a, C 



32y/2A ' 

where we define t = Re T, r = Im T, and treat Wq, A, B as real parameters for simplicity. Note 
that we have dropped the e~ 2x , e~ 2l3x and e~ x ~ l3x terms as well as O ((y) 2 ) terms in anticipation 
of a large x or, equivalently, a large V. This approximation is essentially the same as that used 
for the case with a single non- perturbative term (rj[8j [l2j [l3j [31] . Here, we are interested in the 
parameter region where we can have a large dimensionless volume V ^> 1. Once we have the 
desired solution, we can check a posteriori the validity of this approximation. 



Note that this model (2.5) reduces to the single non-perturbative term model if either (1) 



(3 = 1, where A — > A + B, or (2) B = (z 1 = 0). This single non-perturbative term case has 



been studied carefully 12 (where C = 9/2 C is used) and a meta-stable dS vacuum is obtained 
(with r = 0) only if 

0.811 < C < 0.864, 5/2 < x < 3.11 (2.6) 

where the lower x value corresponds to a Minkowski vacuum and the upper bound indicates the 
vanishing of the modulus mass. This also requires Wq/ A (and so W^o^) to be negative. The 
validity of the large volume approximation requires 

£/2V ~ a 3 / 2 £/12V6 < 1. 



The a'-correction term £ is given as in (2.1), and we do not expect this term to be tiny since 
we like to stay in weak coupling regime Re S > 1 and — \ ~ (9(100). Although we may take 
the volume large by taking a smaller a, we do not expect this to happen naturally. The non- 
perturbative term is obtained by Euclidian D3-brane or gaugino condensation on D7-branes. 
For instance, a = 27r/N for SU(N) gaugino condensation on D7-branes. Recently it is analyzed 



that the D7-brane tadpole cancellation 32 as well as the holomorphicity of D7-branes suggests 



an upper bound for the maximal rank of gauge group on D7-branes 31 , 33] (where N < 24 or 
a > 0.26 in a specific model) . So there exists an upper bound on V itself for a non-racetrack 
type of the Kahler Uplift scenario. 



There are a number of meta-stable solutions to V (2.5): 



For very small B, the B term may be treated perturbatively. Then the single term con- 



straints (2.6) will only be slightly modified. That is, the volume V is still strongly con- 



strained from above. 
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For 3> 1, the dropped e~ 2a; term will become more important than the kept e~ l3x term 
in the approximate V (2.5). In this case, we should restore the e~ 2x term and may treat 
the B term perturbatively. Again, the volume V will be strongly constrained as before. 

For small 0, it turns out that -C 1; in fact, ~ l/sc, and so the dropped e~ 2 ^ x term will 
be more important than the kept e~ x term. Assuming we have a solution when ~ l/x, 
then the volume is constrained above by V ~ (iVfc/27r) 3 / 2 < 21 where we used the bound 
Nb < 24 for Pf 116g 31 as an illustration. 

Since we like to find solutions where x may be taken large without much constraint, we 
shall focus on the solutions where ~ and B needs not be small. 



2.3 Large volume solution 

Since the r ^ case at large volume has no stabilized solution, as explained in appendix [A] 
let us focus on the y = ar = case. Here, we like to show that the stability conditions can 
be met at the large volume approximation for y = and ~ 0(1). Interestingly, the solution 
with > 1 and large volume imply a classically stable dS vacuum with an exponentially small 
cosmological constant A. Also we shall check that the solution of the approximated potential 



(2.5) is well justified even in the full potential (2.3). 



Now we solve for the classically stable minima in the ~ region. The extremal 

conditions dtV = d T V = may be expressed as the relations: 

1 w _ 1)x -2x + 5 + 9^X 7/2 ((3-l) + e x (x 2 + 2x)\ 

z /3(2/3x-5) ' 2/3x-5 ' 1 j 

where the rescaled potential A will be constrained by the stability condition. Plugging z, C in 
the Hessian didjV oc didjX and taking the large volume expansion x ^> 1, we get 

= m 2 > 0, 

_.. , 0-1 5(0 , {90 45 \ (2. 



=- m2+A (5) + "' so - 

while the off-diagonal component d x d y V\ ext = at the y = value. So the stability condition 
(positive mass squared for both x = at and y = ar at the extremum) puts a strong constraint 
on the value of A = — 2V\ cxt /a 3 AW . 



The stability condition (2.8) in the large volume approximation yields (3 > 1 and 



2(0-12/ 5(0 + 2(0-1) / 5(0 + 1) 3 

6 -9^~ I, 1 " ~2^~) " A ~ e "^0^l 1_ ^0^ + 20^ <2 " :)) 
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which takes the form < A m i n < A < A max . So we see that a positive but small A is guaranteed 
together with the large volume V and (3 > 1, implying that a positive cosmological constant 
A = V\ min = (— a 3 A W /2) A = a 3 \A W \X/2 emerges. It is interesting that the stability condition 
imposes no upper bound on x, in contrast to the single term case (2.6). In fact, A m i n — > A max 
as x — > oo, so A approaches an exponentially small positive value at the large volume (x —> oo) 
limit, 



A 



a 3 \AW n 



, 203-1) 
9f3x 



1 - 



5(/3 + l) 
2/3x 



(2.10) 



As we will analyze the probability distribution -P(A) in the next section, a small A is quite 
generic in the Racetrack Kahler Uplift model. 

It is also interesting to find the lower bound on x. To do so, we have to start with the 
exact formulae for <9j<9jA| ex t- With that, we can easily write down the large volume limit of the 
constraints on A, 

_ x (/3-l)(2/3x-5(/3 + l)) 



< A < e 



a .(/3-l)(4/3x 2 -10(/3 + l)x + 35) 



(2.11) 



We see that A r 
bound for x 



9fi 2 x 2 ~ ~ 9x(2/3 2 x 2 - 3f3x - 10) 

— > oo as its denominator factor {2f3 2 x 2 — 3(3x — 10) — > 0. This yields the lower 

3 + ^ 3.11 



As x approaches this bound, the range of allowed A becomes infinite. For x smaller than this 
value, A max < 0. Note that this lower bound matches the upper bound (2.6) as /? — > 1. 

So far we have restricted the parameters in terms of x. Now we estimate the solution for x 
in terms of the parameters in the large volume limit satisfying all the conditions. Inserting the 
leading order of A (2.9) back into the extremal conditions (2.7), we have at large x, 

2(/3- 



1 



(3 3 



C 



9(3 



-e"V/ 2 . 



(2.13) 



If we solve the first equation for x, we get 

1 



x 



ft 



In 



—z 



C 



9/3(/3 - l) 5 /2 



—z 



-1 

p-1 



(2.14) 



Therefore (3 close to one as well as small \z\ contribute to the large volume. To obtain the large 
volume solution, we consider the parameter region so that < — f3~ 3 z < 1. Plugging the large 
volume solution (2.14) to the leading order of the cosmological constant A obtained in (2.9), the 
absolute values of the potential minimum is given by 



A=V\ 



a 3 BW \ 2{(3-l) 2 (3 2 {~(3' 3 z)T- 



m(-/?-V)' 



(2.15) 
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Note that A in the coefficient is rewritten in terms of z. Since the range of A restricted by 
the stability analysis is quite narrow for a given value of z, the formula here gives the right 



magnitude for A at the meta-stable vacua. Next, let us see how the condition for C (2.14) can 

)er x, as is apparent from the formula in (2.3), the 
assuming weakly string coupling Re S = l/g s > 1. 



be met. If we have a non-trivial Euler number x, as is apparent from the formula in (2.3), the 
smallest possible £ may be of order C(10~ 3 ) 



Together with a lower bound for a, a natural requirement to satisfy C in (2.13) is having a small 



W . So A is further suppressed due to the small W in (2.15). To make this point clear, we 



substitute W$ using equation (2.13) and get 



A 



(2.16) 



When P > 1 and simultaneously \z\ is small, we have an exponential suppression in the cosmo- 
logical constant. 

Finally let us estimate how the a'-correction is controllable. Using the large volume solution 



(2.7) and (2.14), we get 



C 



(-/r 3 *)^- 1 |m(-/rV)| 



(2.17) 



x 3 / 2 9(3((3 - 1) 

Since this uplifting term is highly suppressed as a function of z and /3, our approximation 



keeping up to the linear term for a'-correction in (2.5) works quite well. On the other hand, the 
suppression ratio for the a'-correction becomes 



J_ a 3 /^ (/3-1) 3 / 2 
2V ~ 4^2 |ln(-/3" 3 z)| 3 / 2 ' 



(2-18) 



It is clear that if we have only single non-perturbative term, there is no suppression depending on 
z, (3. So this large volume approximation also works to make the required a'-correction smaller 
such that type IIB SUGRA approximation stays valid. The construction with the large volume 
makes several corrections under control, including higher a' or stringy loop corrections, which 



may scale as 0(V 10/ ' 3 ) in the potential 15-20 in light of the extended no-scale structure 17,18 



Here we present an example for illustration. Let us start with the following input parameters: 



2tt 
15 : 



Wo = -0.223, A = 1.65, B = -4.77, a 
So the combined parameters are given by (to 3 digits) 
z = -0.346, C = 8.28 x 10~ 6 , (3 = 



15 
14 



2tt 
14' 



1.07, 



£ = 3.41 x 10 



-3 



0.281, 



(2.19) 



(2.20) 



We find the minimum at t — 42.5, r = in the full potential (2.3) as shown in figure [TJ where 
the volume is quite large, V = 785. 
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42.8 



Figure 1: The minimum in the full potential (2.3) with the parameters given in (2.19). The minimum sits at 
* = 42.5, r = 0. 



The solution in the full-potential (2.3) with input (2.19) is given by 



x 



(Mi) = 17 _ 82 28, A = U (full) | 



1.60122 x 10 



-13 



while the approximate potential (2.5) yields 



x 



(approx) = 17 _ 82 1Q y (^Pprox) | = L60 H 5 x 1Q- 13 . 



On the other hand, the approximated analytical formulae (2.14) and (2.16) yield 

x (anal) = ^7597 y(anal), = L78635 x 1Q - 



-13 



(2.21) 



(2.22) 



(2.23) 



So we see that the approximations work quite nicely not only for the approximate potential 



(2.5), but also for the analytic expressions. 



3 Probability distribution of Racetrack Kahler Uplift 



In previous section, we see that the racetrack Kahler Uplift model have no upper bound for the 
volume moduli, unlike the Kahler Uplift model with a single non-perturbative term. To under- 
stand how likely a tiny cosmological constant will appear in this racetrack model, we analyze 
the probability distribution -P(A) in this section. The A at the classically stable minimum of 



the potential (2.16) is given by 



_ 62y/2a 3 / 2 i? 2 /3 4 « 5/2 

A = kp- 1 — hiK 

243£v/^T V ; 



K = 



(3.1) 
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Now we introduce a randomness to the system. As discussed in [7}]9| [34|[35] , when we take 
into account the moduli stabilization of complex structure moduli and dilaton with different 
values of fluxes, we expect that the different values of Wq, A, B are given. Together with the 
fact that there are also many types of models for complex structure moduli, corresponding 
to many varieties of Calabi-Yau compactifications, we have a rich enough structure of vacua 
in the string landscape. To deal with all of these models is rather complicated, so we mimic 
this variety by just simply randomizing some parameters. Let Wq be a random parameter 
obeying a uniform distribution with a range that satisfies the condition for C in (2.9), so the 
expression (3.1) follows. To simplify the analysis, we set B — 1 and £ = 1. Since B contributes 
just in the coefficient and does not really touch the details of the dynamics, we set simply 
B = 1 for simplicity of the arguments. Learning from the analysis in (7j[8], it is clear that 
randomizing B will probably not diminish the divergent peak in the distribution P(A) towards 
A = 0. We also do not randomize the parameters a, (3. The Euclidian D3-brane gives us the 
non-perturbative term with a = 2tt, while the gaugino condensation on D7-branes produce the 
term, e.g. with a = 2tt/N for SU(N) group. But since we also have to satisfy the tadpole 
cancellations in Calabi-Yau compactification, which affects the number of the D7-branes to keep 
its holomorphicity 31 33 . So since there may not remain large choices for these a, b, we rather 
pick up a value for (5 and set a = 1 for simplicity. 

Note that the statistical distribution of the flux vacua was considered by giving a distributed 
randomness to the flux quantities in 36-38 . In contrast, our interest here is to estimate the 
probability distribution P(A) of A of the flux vacua, especially in the concrete dS cases so that 
the stabilization dynamics crucially affects the distribution. The shape of the distribution would 
imply how likely we can achieve small values for A. 

By setting a = B = £ = 1 for simplicity, the probability distribution is estimated by the 
formula: 



62^2/3 



4 



P(A) = j dz P{z) 5 \A - (- ^ ^1 ■ (3.2) 

For a fair discussion, we consider the uniformly distributed — 1 < z < with P(z) = 1 as a 
conservative choice. Then the integration can be performed quite easily by 

243Q? -!)»/» *n 

1 ] 32x/2/3 (-lnK) 3 / 2 (5-5/3-4/31nK)' K ] 

The range — 1 < z < is good for the large volume approximation. For instance at /3 = 1.1, 
we have x ~ — ln«/(/3 — 1) > 2.86, and the distribution (3.3) is a well-defined monotonically 
decreasing function of k. 



Since we would like to rewrite this as a function of A, we solve (3.1) for k by 

'-=^ w -'(- 5X 23/X-ir/j - (3 ' 4) 
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P(A) 



Pfln A) 



fixlO 1 

5 xlO 1 
4X10 1 
3 xlO 1 
2X10 1 
1 xlO 1 



0.04 
0.03 
0.02 
0.01 



. xlO" 



4. xlO" 11 4. xlO" 11 I. xlO" 



1.x 10"^ 




-151 



-100 



-iO 



In A 



Figure 2: The asymptotics at ft = 1.1 of the analytical expression (3.5| and the numerical data. 



Here we introduce the Lambert VK-function which is the solution ofWe w = X, and W = W-i(X) 
when W < — 1, since our large volume approximation works if 4/3 In k/5(/3 — 1) ~ — 4/3x/5 <C — 1. 



Inserting this solution into (3.3), we get 



P(A) 



W_ 1 =w_i 



loovlocs - 1) (-W_ 1 )3/2(l + w_i) ' (3 5) 

9A 2/5 



5 x 2 3 / 5 /3 3 / 5 (/3 - I) 4 / 5 J ' 



Let us expand (3.5) around A ~ to get the asymptotic behavior. Using the expansion of 
W_i for small X: 

W-i (X) ~ In X - ln(- In X) + ■ ■ ■ , (3.6) 
the probability distribution becomes 

16(0- 1) A^(-lnA) 5 / 2 

So for /3 > 1, we see that the diverging behavior is very peaked as A — > 0. Since (/3 + 1) / 2/5 < 1, 
P(A) is normalizable, i.e.,/ P(A)dA = 1. 



We illustrate the result in the figures of the probability distribution function of A in (3.2). 



Here we again set a = B = £ = 1 and uniformly distributed — 1 < z < 0. The analytical 



expression (3.5) as well as the numerical histogram of (3.1) are illustrated in figure 2 at = 1.1. 



The distribution is quite sharply peaked toward A = 0, as estimated in (3.7). 

To get a better feeling of quantification of the peaking behavior, it would be better to 

/o 



introduce the likely value that Y% of the data fall in the value Ay: f Q Y <iAP(A) = Y%. Using 



the data obtained above at j3 = 1.1, the likely values become 

A10 = 3.61 x lCT 24 , A 50 = 7.08 x 1(T 10 , A 80 = 4.25 x lCT 6 , (A) = 8.90 x 10~ 6 , (3.8) 
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where we also present the average value for comparison. We see that just 10% fine-tuning 
suggests the substantial suppression of cosmological constant. This is nothing but because 
the highly sharply peaked behavior as in (3.7). Note that the A 50 is simply the median. For 
comparison, we have, at /3 — 1.04, a much sharper peaking behavior as is clear from (3.7): 

-19 a/9=1-04 



A 



/3=1.04 
10 



2.83 x 10 



-54 



» ,3=1.04 
iv 50 



5.47 x 10 



NO 



2.81 x 10" 



=1.04 



6.36 x 10" 



(3.9) 



4 Swiss-Cheese type model 

So far we have focused on a single Kahler modulus case. Here we introduces multi-Kahler moduli 
and check whether the multi-moduli case is compatible with the large volume approximation in 
the Racetrack Kahler Uplift, especially with the Swiss- Cheese type of compactification. 

The Swiss-Cheese type model is a class of Calabi-Yau compactification, and is used to realize 
the Large Volume Scenario (LVS) [14) . It is clarified that there is a large variety of Swiss-Cheese 



type of compactification 19 39 40 . In LVS, the volume is made actually quite huge. The 



large volume is good to have a control of several corrections including higher a' or stringy loop 



corrections, scaling like 0{V~ 10 ^ 3 ) in the potential 15 20 . In our analysis with the Racetrack 
Kahler Uplift, although we do not expect that the volume can be huge naturally as well as LVS, 
we consider that the volume may be large enough so the corrections are under control. 

We focus on a two Kahler moduli test example to investigate the multi-Kahler 

scenario. The model is given by 

K = -2ki fv+|) +■■■ , 

V = (T 1 +T 1 ) 3/2 -(T 2 + T 2 ) 3/2 , W 

W =W + A ie ~ aiTl + B ie - blTl + A 2 e- a2T \ 

Here we just introduce single non-perturbative term for the second modulus. Again we assume 
that the complex structure moduli and dilaton are stabilized supersymmetrically, and choose 
the solution for imaginary modes to be Im T\ = Im T 2 = for simplicity. We are interested 
in the parameter region which include dS vacua as a result of the precise Kahler uplift, the 
potential may be approximated up to leading order of the non-perturbative term as well as the 
a'-correction, by 

V =e K (\DW\ 2 -3\W\ 2 ) 

a^B 1 W z\ ( x\e~ xi /3xie _/3xi , nS , x 2 e~ X2 C 

cosj/i — cos(/fyi) - d 2 — cosy 2 1 



2 J I v 2 z v 2 " v 2 v 3 / ' (4-2) 

\ 3/2 A 

x 2 \ A 2 a 2 

— , x 2 = a 2 t 2 , y 2 = a 2 r 2 , d 2 = —, d 2 = —. 
b 2 A 1 a 
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We analyze the dynamics of this effective potential. 

First, we consider the dynamics on y 2 direction, where the second derivative is given by 

2 d 2 x 2 e~ X2 
d y2 V=— 2 cos ?/ 2 . (4.3) 

Together with the fact that the first derivative is proportional to sin 2/2, we can have a stable 
solution at ?/2 = when d 2 > 0. The off-diagonal component with respect to y 2 at the extrema 
didy 2 V\ ext is now trivial due to the solution y 2 = 0. Since yi = solution is also motivated as 
argued in the previous section and x 2 direction can not change the dynamics for xi,y% at the 
large volume, so we take the solution with 2/1 = 2/2 = 0. 

Similarly to the previous analysis, the extremal condition d^V = with A = —2V\ ext /a 3 BW z 
can be rewritten by 



z \ P 2 2/3%! J + V2/3 2 + 4/33 

e -xi+x a _ i5(/3_ i)a;V 2 \ e ^ A / 9^/2^1/2 2 7x\ /2 x 1 2 /2 ' 

2 ~ 1 - x 2 y 2l35 3 2 /2 x{ /2 Af5 2 5l ,2 x^ 2 J l-x 2 \ 25 3 2 /2 A/35 3 2 /2 



Plugging the extremal solution into the Hessian, we get, at the large volume limit 
afB.WozY 1 o2 „. 5(/9 — 1)\ x / 9/3 9 



(4.4) 



2 , d ^^~ x \^-^) + \-^-M ] + 

a\Br W 0ZV > ^ 5(/3-l) X + A f 9/3 + 45 



2 ; yi ext V ^? 2x? ; T 

alBxWpz y 1 2 _ gi 3(/3 -!)(! + 3x 2 -2x 2 ) 9(1 + 3x 2 - 2x|) 

O J ^a 17 lext ~e .^,3/2 7/2 1/2/, n + A , <-3/2 3/2 1/2/., 

z / 4/302' x{ x 2 (1 — x 2 ) 4o 2 x{ x 2 (1 — x 2 



a\B x W,z \ x 9(/3-l)xf 81xf 

2 J <V^U~ 2 /35fxf 4/35 2 3/ M /2 ' 

(4.5) 

where we keep the next-leading order terms in d 2 x ^/\ cx t,d 2 u ^V\ e ^\, since their leading order terms 
may vanish due to the stability condition, while the leading order terms in d% 2 V\ ext , c^i^^lext 
are not. 

According to the Sylvester's criteria, the positivity of the sub-matrices are necessary con- 
ditions for the positivity of the entire matrix (see e.g. El], also applied for necessary stability 
constraints in 42-44]). So we consider the stability in the X\-y\ subspace first. Similarly to the 
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previous section, the condition in this subspace may be expressed as 

\ 9pxi 9p z x\ J \ 9px\ 9p 2 xf J 

for large x\ and (3 > 1. Recall that d 2 > is necessary for the stability on the y 2 -direction, 
which becomes 



0<.'- 2 <l. -e~ Xl ^4<\, or l<x 2 , A<-e-^f-4. (4.7) 
3px( 3px( 



So, to meet the condition (4.6), the solution is required to stay within < x 2 < 1. 

The remaining task is to check the stability condition in the x\-x 2 subspace. Now plugging 
the leading order of A into the Hessian, we get 

alB.Woz y 1 2 - 8l 3(/3-l) 
2 J 9 ^~ e - 1 ^- r + -, 

a^Wpz y 1 a2jr[ _ xi (/3-l)(l + 3x 2 -2x^ 

2/35fxfxf (l-x 2 ) 
^■nn-. y- 9(/3 2 -l)xf 



2 y " + -» ^ 



It is clear that the diagonal components are positive for > 1, < x 2 < 1, while the off-diagonal 
components are sub-leading in the determinant; therefore, positivity of the Hessian is assured 
at the large volume limit. 

Let us summarize the stability analysis in this section. The extremal conditions for meta- 
stable vacua become 

-"—eV- 1 *", d 2 ~e-* ^ 6~e^ 2 ^~ l } X *\ P > L (4-9) 

z /3 3 p5l /2 {l-x 2 ) 9/3 

Then the solutions sit in the region 

1 < xi, < x 2 < 1. (4.10) 

So having the racetrack type of potential for big volume modulus X\ is well-motivated to realize 
the large volume even at dS vacua. Since we can easily have large volume solutions at dS vacua, 
we can control the several stringy corrections, simultaneously realizing the cosmological constant 
which scales exponentially small: 

A^|^^(- r ^(-M- r - 2 ,f, (4.11) 



where we have used the approximate solution for x\ using (4.9), similar to (2.16). 
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2.10 



Figure 3: A dS minimum in the full potential (4.1) 



We end this section by presenting a sample dS solution. In figure |3j we illustrate a dS 



solution in the full-potential (4.1) with the parameter set: 



W = -0.223, Ax = 1.65, B x = -4.77, A 2 = 2.24 x 1(T 7 , 

2tt 2tt 2tt , ( 4 - 12 ) 

ax = — , bx = — , a 2 = — , £ = 4.35 x 1(T 3 . 
15 14 15 

Then the solution is given by 

ti ~ 42.5, t 2 ~ 2.18, n = r 2 = 0. (4.13) 
So the large volume is easily realized: V ~ 774. At the local minimum of the potential, 

A = 1.21 x 1(T 13 (4.14) 

which is close to the previous result. 
The Hessian is given as 

dlV\ min =1.07 x 10- 15 , dlV\ min = 2.66 x 10" 14 , d tl d t ,V\ mi n = -1.33 x lO" 16 , 

(4.15) 

d 2 T1 V\ min =3.79 x 10- 16 , d 2 T2 V\ min = 2.14 x 10" 14 , d T1 d T2 V\ min = 1.17 x lO" 20 . 

Therefore the minimum is stable in general. There is no other cross term between the real and 
imaginary parts owing to n = r 2 = 0. 

Note that the combined parameters here are given by 

z = -0.346, d 2 = 1.36 x 10" 7 , C = 1.06 x 10~ 5 , (3 ~ 1.07, 5 2 = 1., 



(4.16) 

x x = 17.8, x 2 = 0.915, y 1 = y 2 = 0. 
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5 Discussions 



In this paper, the possibility of the large volume X\ 1 is allowed owing to the racetrack 
potential. Since the volume is a parameter to control the a'-corrections and also the string- 
loop corrections of order (9(V~ 10 / 3 ) in the potential 15-20 in light of the extended no-scale 



structure 17, 18 , the realization of the large volume is well-motivated to achieve a dS vacuum 
which is meta-stable even in the presence of corrections. The previous Kahler Uplift model is 
basically constrained by the dynamics, and suggests x\ ~ 3. Even though the large rank of 
the gauge group on D7-branes helps to relax the upper bound on the volume, we should be 
concerned with this way of relaxation due to the constraint of the gauge group rank by the 



D7-tadpole cancellation condition 31,33 . So relaxing the constraint for x\ in the racetrack 
model is important for the construction of dS vacua. 

It is worth mentioning that the resultant cosmological constant is exponentially suppressed 
as a function of the volume. At a dS vacuum with a small cosmological constant, the corrections 
would be more suppressed owing to the large volume. However, we should keep in mind that the 
combined parameter C is required to be exponentially small as well. C can be small due to an 
exponentially small W . In fact, the peaked distribution of W toward W = is obtained using 
a linear model for complex structure moduli |8]. We may expect that the sharper peaked distri- 
bution of Wo is realized in the presence of more non-trivial couplings. So the even exponentially 
small Wo is quite conceivable. 

Recently, a new a'-correction is estimated for M = 1 compactification when the first Chern 
number of a three-dimensional Kahler base B 3 in M/F theory is non-trivial 45 . Since the coef- 
ficient for this correction appears non-negligible, the volume is required to be large to suppress 



this correction 46 . Here, the important suppression parameter for this additional correction 
would be the ratio of the coefficients between y~ 10 / 3 term owing to the extended-no-scale struc- 
ture 17,18 and the leading a' correction term [25] scaling like V -3 in the potential. Including 



this r\ correction term to the potential V (2.5) we study in this paper, we would have (V ~ x 3 / 2 ), 

J, J_ _ V_ _ tg s =i _ V_ 

x 9/2 x 9/2 x 5 x 9/2^/2 x 5 ' 



where rj is the coefficient of the new a'-correction and the string coupling dependence of £ (2.1) 
is made explicit. Our qualitative result will remain valid if this remains positive for relatively 
large volume. This can be satisfied if either the volume is large or if rj < £. This later condition 
may be satisfied in the weak coupling approximation. Of course, the Racetrack Kahler Uplift 
model analyzed in this paper is applicable for large classes of models including those with trivial 
first Chern number of P3, in which case rj = 0. 

The probability distribution P(A) is sharply peaked toward A = 0, explaining a natural 
statistical preference for a small cosmological constant. The distribution P(A) is diverging, but 
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normalizable. Since this mechanism works to have a hierarchical structure from Planck scale, we 



may worry about the cosmological moduli problem 47 49 . The cosmological moduli problem 
is a constraint from reheating of the universe, so it requires details of how the moduli fields 
decay to the matters. Also there are some ways to relax the constraint, including a thermal 
inflation which is a mechanism to dilute the energy produced by the moduli coherent oscillation. 
Recently, it is analyzed in detail that a double thermal inflation relaxes the constraint in the 
Large Volume Scenario 50 . Therefore, the cosmological moduli issue crucially depends on the 
detail of cosmological history of the universe. 
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A Case for non-zero axion value of the Kahler modulus 



Here we like to show that the stability condition of the system with y ^ at large volume cannot 
be satisfied. We first analyze the case for non-trivial imaginary mode y ^ in the approximated 



potential (2.5). Since d y V = is automatic for y = 0, we shall assume y ^ 0. The extremal 
condition d x V = d y V = give the relations 

x 5/2 

C = e~ x — {(3{x + 2) cosy - {(3x + 2) sinycot(^y)) . (A.l) 



1 

z 



,08-i)x sm v 



/? 2 sin(/?y)' ~ 
Substituting z, C in the potential, we get 



V 



ext. 



a 3 B W z 



x>1 ( a 3 B W z 



9(3x 2 
, 2 



(Ji(2x - 5) cosy + (-2#c + 5)siny cot(/3y)) 
cosy — siny cot (/?£/)) + ■ ■ ■ 



(A.2) 



2 J 9x 

where in the last equation, we took the large volume approximation assuming j3 ~ 0(1). 



Next we consider the classical stability condition. Again, using the relations ( A.l ) and taking 
the large volume approximation for V (2.5), components of the Hessian at the extremal didjV\ ex t 
are approximated by 



d 2 x V\ ext 



a 3 BW z\ _„ 1 

e — (cosy — smy cot (pyj) + 



X dyV\ 



ext 



a 3 B W n z 



a 3 B W n z 



(cosy — siny cot (/3y)) 



(A.3) 



,0 



x- 



smy. 
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Since the sign of the leading term of c^V| ext is opposite to dyV\ ext , we see that the leading term 
of the determinant is negative, i.e., det (didjV\ ext ) < 0. Note that the form of dyV\ ext is exact 
without the large volume approximation. So no matter which cosmological constant we have, 
we cannot satisfy the stability condition for the system with y ^ at the large volume. 
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